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Eddington ratio governs the equivalent width of Mg n emission line in active 

galactic nuclei 
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ABSTRACT 

We have investigated the ensemble regularities of the equivalent widths (EWs) 
of Mgll A2800 emission line of active galactic nuclei (AGNs), using a uniformly 
selected sample of 2092 Seyfert 1 galaxies and quasars at 0.45 < z < 0.8 in the 
spectroscopic data set of Sloan Digital Sky Survey Fourth Data Release. We find a 
strong correlation between the EW of Mgll and the AGN Eddington ratio (L/L Edd ): 
EW(Mgll) oc (L/LEdd) -0 ' 4 • Furthermore, for AGNs with the same L/L Edd , their EWs of 
Mg II show no correlation with luminosity, black hole mass or line width, and the Mg II 
line luminosity is proportional to continuum luminosity, as expected by photoioniza- 
tion theory. Our result shows that Mg II EW is not dependent on luminosity, but is 
solely governed by L/L Edd . 

Subject headings: radiation, accretion - galaxies: active - quasars: emission lines - 
quasars: general 
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1. Introduction 

Mg II A2800 is one of the most prominent broad emission lines in the near-ultraviolet spectra 
of type-1 active galactic nuclei (AGNs, including Seyfert galaxies and quasars). Its equivalent 
width (EW) is generally in the range of about 10-100A (Dietrich et al. 2002). Theory shows that 
as a low-ionization line, Mgll originates from optically thick (i.e. large column density) clouds 
only that suffer little from radial motions, thus it is valuable in probing the properties of the clouds 
gravitationally bound in the AGN broad-line region (BLR), and in estimating the virial mass of the 
central supermassive black holes (e.g. Grandi & Phillips 1979, McLure & Dunlop 2004, Wang et 
al. 2009). 

However, Mg II is highly blended with the Fe II multiplet emission. This complicates accurate 
Mg II measurement in AGN spectra, and makes a comprehensive study difficult of Mg II using 
large, homogeneous AGN samples. By composite-spectrum analysis, Dietrich et al. (2002) found 
a negative correlation between the Mg II EW and the continuum luminosity, similar to the Baldwin 
(1977) effect first discovered in ClV line. However, except that and some others using small 
samples (e.g., Grandi & Phillips 1979, Zheng et al. 1992, Espey & Andreadis 1999, Croom et al. 
2002), there has been few studies so far on the Mgll properties using a large, homogeneous AGN 
sample. 

Recently, this technical difficulty has been overcome by the availability of the detailed Fe II 
templates in the wavelength region covered by Mgll (e.g., Tsuzuki et al. 2006; cf. Vestergaard & 
Wilkes 2001). In this Letter, we explore the regularities of the Mgll EW in the AGN ensemble, 
by taking advantage of the unprecedented spectroscopic data from the Sloan Digital Sky Survey 
(SDSS; York et al. 2000) and the new Fe II templates. We find that, at the Oth-order, Mg II luminos- 
ity is directly proportional to continuum luminosity; at the lst-order, the proportional coefficient 
correlates negatively with the Eddington ratio (£ = L/L E dd)Q Throughout this paper, we use a 
cosmology with H =70 km s" 1 Mpc" 1 , Q M =0.3 and il\=0J. 



'Eddington ratio (I = L/L^dd) is the ratio between the bolometric and Eddington luminosities. Eddington lumi- 
nosity (LEdd)> by definition, is the luminosity at which the gravity of the central source acting on an electron-proton 
pair (i.e. fully ionized gas) is balanced by the radiation pressure due to electron Thomson scattering. Thus Eddington 
ratio is different from relative/normalized accretion rate both in meaning and in scope of application; see Dong et al. 
(2009) for a detailed discussion. For clarification and for the ease of notation, below we use L/Lem or I alternately. 
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2. Sample and Data analysis 

We use the UV sample of 2092 Seyfert 1 galaxies and quasars of Dong et al. (2009), selected 
from the spectral data set of the SDSS Fourth Data Release (Adelman-McCarthy et al. 2006). 
Sample definition and data analysis methods are described in detail in that work. Briefly, we 
select quasars with both the H/3 and Mgll present in the SDSS spectra and with continuum and 
emission lines suffering little from the contamination of the host-galaxy starlight. The criteria are: 
(a) 0.45 < z < 0.8, (b) the median signal-to-noise ratio (S/N) > 10 per pixel, and (c) the rest- 
frame absorption-line EWs of CaK (3934A), CaH + He (3970A) and H5 (4102A) < 2a. The H/3 
BLR has been extensively studied using reverberation mapping for about forty type-1 AGNs, and 
the black hole mass formalisms based on single-epoch H/3 are well calibrated with reverberation 
mapping data (Peterson et al. 2004, Bentz et al. 2006, Bentz et al. 2009). Thus the presence of H/3 
in our sample is helpful to double check the findings in this Letter (see §4). 

As described in Dong et al. (2009), the SDSS spectra are corrected for Galactic extinction and 
the Mg II emission is separated from AGN continuum, Balmer continuum and the Fe II multiplets. 
Then the broad components of the Mgll AA2796, 2803 doublet lines are each modeled with a 
truncated 5-parameter Gauss-Hermite series profile (see also Salvainder et al. 2007). The narrow 
component of each line is fitted with a Gaussian. In the fitting, the narrow component is constrained 
as follows: FWHM < 900 km s" 1 and flux < 10 per cent of the total Mg II flux (Wang et al. 2009; 
cf. Wills et al. 1993). UV Fell is modeled with the tabulated semiempirical template generated 
by Tsuzuki et al. (2006) based on their measurements of IZW 1; in the wavelength region covered 
by the Mgll emission, this template uses the calculation with the CLOUDY photoionization code 
(Ferland et al. 1998). 

All the data used in this Letter are presented in Dong et al. (2009). The detailed fitting results 
are available online for the decomposed spectral components (continuum, Fe II multiplets and other 
emission lines) j_ We compared the Mg II parameter values fitted by our fitting method with those 
by the methods of McLure & Dunlop (2004) and Salviander et al. (2007). Our results are well 
consistent with those obtained by Salviander et al., and roughly consistent with those obtained by 
McLure & Dunlop (Wang et al. 2009; cf. Salviander et al. 2007). 

We estimate the measurement uncertainties of the parameters using a bootstrap method (Wang 
et al. 2009; cf. Dong et al. 2008). The estimated typical la relative errors are 10% and 8%, 
respectively, for the fluxes of broad Mgll and H/3; 20% and 15%, respectively, for the FWHM 
of broad Mgll and H/3; 8% and 5%, respectively, for the slope and normalization of the local 



2 Available at http://staff.ustc.edu.cn/~xbdong/Data_Release/ell_effect/, together with auxiliary code to explain and 
demonstrate the fitting and the parameters. 
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continua. According to standard propagation of errors, the formal la error of Mgll EW (in log- 
scale) is typically 0.05 dex. We must note that the thus estimated errors still do not account for the 
uncertainties caused by Fe II template mismatch, etc., and that it is almost impossible to pin down 
the real measurement errors for individual objects. Hence we waive the analysis on the intrinsic 
scatter of the relations of interest in this Letter. 

We calculate the black hole masses, M BH , based on the Mg II FWHM and the monochromatic 
luminosity L 30 oo = AL A (3000A) using the formalism of McLure & Dunlop (2004). The small dis- 
crepancy of the fitted Mg II parameters between our method and the method of McLure & Dunlop 
(2004) has no significant effect on our conclusions (see also §4). The Eddington ratios are calcu- 
lated assuming a bolometric luminosity correction L bo \ fa 5.9L 30 oo (McLure & Dunlop 2004). The 
mean and standard deviation (computed in log-space) of Mgll FWHM are 3200 km s" 1 and 0.18 
dex; L 30 oo, 1.2 x 10 45 ergs s" 1 and 0.31 dex; M BH , 1.6 x 1O 8 M and 0.42 dex; L/L Edd , 0.39 and 0.38 
dex. 



3. Results 

3.1. Correlations of Mg II EW with other quantities 

We use the Spearman rank method to perform bivariate correlation analysis of Mg II EW with 
Mg II FWHM, L3000, M B h and L/L Edd . The results are summarized in Table 1 (top panel), giving the 
Spearman coefficient (r s ) and the two-tailed probability (P nu u) that a correlation is not present. The 
strongest correlation of Mg II EW is with L/L Edd (r s = 0.59) which is a combination of Mg II FWHM 
and L3000 as L/L EM oc FWHM" 2 !'^,,. The correlation of Mg II EW with another combination, M BH 
(oc FWHM 2 ^,), is much weaker (r s = 0.38). The correlation with FWHM (r s = 0.55) is slightly 
weaker than that with L/L Edd , while the correlation with L 30 oo is weak (r s = 0.2). These correlations 
are illustrated in Fig. 1. 

We must note that because the SDSS spectroscopical survey is magnitude-limited, Mgll 
FWHM, L3000, ^bh and L/L Edd in our sample correlate one with another apparently. The apparent 
(likely not intrinsic) correlation between M BH and L/L Edd is further enhanced by the correlation of 
their measurement errors, because both M BH and L/L Edd are constructed from Mgll FWHM and 
L 30 oo- The Spearman correlation coefficients (r s ) of L/L Edd with Mgll FWHM, L 30 oo and M BH are 
-0.92, 0.27 and -0.64, respectively. 

In light of the serious inter-dependence among these four quantities, the correlations of Mg II 
EW with the first three quantities is probably only the secondary effect of the strongest (thus pre- 
sumably intrinsic) correlation with L/L Edd . To test this possibility, we perform partial Spearman 
rank correlation analysis (Kendall & Stuart 1979). The results are summarized in Table 1. When 
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controlling for L/L Edd , none of the correlations of Mg II EW with FWHM, L 30 oo or M Bll is signif- 
icant at P null < 0.01; yet controlling for these three quantities, the correlation of Mgll EW with 
L/L Edd is still significant with P null < 10" 38 . This fact suggests that L/L Edd is the sole parameter, at 
least among the four, that governs the EW of Mg II. 

From Fig. 1 there appears a linear relation between Mgll EW and L/L Edd in log-log scale. 
We perform linear regression, treating L/L Edd as the independent variable, and find: 

logEW(Mgll) = (1.26±0.02) + (-0.39 ± 0.03) \og£ . (1) 

The errors on the fit parameters are purely statistical, and the actual uncertainty could be somewhat 
larger (cf. §2). 



3.2. The dependence of Mg II EW: not on luminosity, but on Eddington ratio only 

As suggested by the above partial correlation tests, the Mgll EW likely does not depend on 
the continuum luminosity or black hole mass, but on L/L Edd only as in Eqn. 1. That is, (1) the 
emission-line luminosity, L Mg n, is directly proportional to the continuum luminosity at a fixed 
L/L Edd , while (2) the proportional coefficient is different for objects with different L/L Edd . 

To further investigate the Mg II EW-L relation, we take an alternative approach to analyze this 
problem. In Fig. 2 we plot the distribution of L MgII versus L 300 o using the objects in our sample, 
with the 433 objects of L/L Edd > 0.8 (high-£ subsample) denoted in blue and the 456 objects of 
L/L Edd < 0.2 (\ow-£ subsample) in red. It is evident that the objects of the high-£ subsample are 
systematically smaller than those of the \ow-£ subsample in L MgII at any given L 30 oo- We perform 
linear regression in log scale with L 30 oo being taken as the independent variable, according to 

logL Mg n = HogL 30 oo + b . (2) 

We fit the relation for the full sample, the high-£ and low-£ subsamples, and the intermediate-^ 
subsample with 0.2 < L/L Edd < 0.8, respectively, with and without k being fixed to be 1. The 
results are summarized in Table 2. The \ow-£ subsample has 456 objects clustering mostly in a 
narrow range of luminosity, yet with a fairly large dynamical range in L/L Edd of about 1 order of 
magnitude. The fit with k being an additional free parameter for this subsample is not statistically 
guaranteed according to F-test. Thus in Table 2 we list only the fitting result with k fixed to be 1 . 
Instead, we select the subsample of 299 objects with 0.1 < L/L Edd < 0.2 and perform the fitting. 
The results are listed in Table 2 too. 

The best-fit k values for the high-£ and intermediate-^ subsamples are well consistent with 1 
within lcr. For the 0.1 <L/L Edd < 0.2 subsample, the best- fit k is consistent with 1 within 1.5a. As 
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Mgll EW is equivalent to L Mg n/^30oo 5 such a direct proportional relation between L Mg n and L 30 oo 
means that the Mg II EW does not depend on luminosity. 

The mean values of logL/L Edd of the full sample, the high-£, intermediate-^, 0.1 < i < 0.2 
and \ow-£ subsamples are -0.41, 0.08, -0.38, -0.83 and -0.95, respectively. The coefficients (b) 
of these best- fit proportion (k = 1) relations are consistent with the predictions of Eqn. 1. 

The best-fit k (0.91 ±0.01) for the full sample is significantly lower than 1, which means 
an apparent dependence of Mgll EW on luminosity for the full sample, as illustrated in Fig.l. 
However, this can be naturally explained by the selection effect mentioned in §3. 1 : high-£ (and thus 
small-EW) objects dominate the high-luminosity end while \ow-£ (and thus large-EW) objects tend 
to cluster in the low-luminosity end, which is clearly displayed in Fig. 2. As a demonstration (see 
also §4), in the panel b of Fig. 1 we over-plot the inferred Mgll EW-L 30 oo relation, the Baldwin 
effect, that corresponds to the fitted logL MgII -A:logL3ooo relation where k = 0.91. 



4. Discussions and conclusions 

To check possible effects of black hole estimation on the above correlations, we re-examine 
the above correlation tests by calculating M B h using all the other available formalisms based on 
broad H(3 and Mg II lines (see a compilation of the formalisms in McGill et al. [2008] and the 
formalisms in Wang et al. [2009] and Vestergaard & Osmer [2009]). All the tests give similar 
results to those listed in Table 1 (top panel). This is mainly because the M B h dynamical range 
covered by our sample is not large 1.5dex), and the various formalisms based on single-epoch 
HP or Mg II have only subtle differences one from another (Wang et al. 2009). In the bottom panel 

o 

of Table 1, we also list the results using the FWHM of broad H/3, the luminosity AL A (5100A), the 
black hole masses calculated with the formalism presented in Collin et al. (2006, their Eqn. 7) and 

o 

the Eddington ratios assuming L boi « 9AL A (5100A). 

Considering the fairly narrow dynamical range of L/L Edd (la = 0.38 dex in log-scale) and 
the measurement errors, such a significant correlation between Mg II EW and L/L Edd is surprising, 
suggesting a rather tight intrinsic relationship between them. More importantly, Mg II EW is likely 
to be intrinsically solely dependent on L/L Edd but not on continuum luminosity. The dynamical 
range of the Mg II EW of the objects in our sample is similar to that of the composite spectra of 
Dietrich et al. (2002). They reported a Baldwin effect of Mgll as logEW oc (-0.09 ±0.01) logL. 
This slope is exactly consistent with our k value (0.91 ; see Table 2 and Fig. lb) for the full sample, 
which is likely due to the selection effect discussed above. However, the range of our L 30 oo is 
2 x 10 44 -2 x 10 46 ergs s" 1 whereas that of Dietrich et al. (2002) is 2 x 10 42 -2 x 10 47 ergs s" 1 
(converted from AL A (1450A) by assuming fx oc A" 15 ). Our sample lacks objects of relatively low 
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luminosity. Thus, if our above finding holds for the AGN ensemble in the full luminosity range as 
in Dietrich et al. (2002), then the classical Baldwin effect of Mg II is purely a secondary effect of 
the EW-£ relation in combination with the selection effect inherent in any flux-limited sample (see 
the analysis in the last paragraph of §3.2). 

It is intriguing to understand the EW-£ relation as expressed in Eqn. 1 (and Table 2; cf. Baskin 
& Laor 2004, Bachev et al. 2008). First of all, the fact that the Mg II luminosity is directly propor- 
tional to continuum luminosity is just as expected from the photoionization theory. Second, in the 
photoionization picture (Osterbrock & Ferland 2006), the emission-line EW is dependent mainly 
on two parameters, the shape of the AGN continuum and the covering factor of the line-emitting 
clouds, thus the fact that the proportional coefficient is different for different L/L Edd suggests that 
either the continuum shape or the covering factor or both is dependent on L/L Edd . 

Zheng & Malkan (1993) and others have proposed a model in which systematic variation in 
the ionizing continuum shape with luminosity explains the Baldwin Effect, via the softening of 
the accretion disk spectrum with increasing black hole mass. More recently, Korista et al. (1998) 
generated grids of photoionization models to demonstrate this relationship quantitatively, and ad- 
ditionally invoked variations in gas metallicity with luminosity to explain the peculiar weakness 
of the Baldwin Effect in N V A1240. The phenomenological model of Korista et al. relating ion- 
izing continuum shape and gas metallicity with quasar luminosity (or M B h) generally explains the 
Baldwin Effect in almost all emission lines (cf. Dietrich et al. 2002). The remaining question is 
the underlying physical link between the shape of the ionizing continuum and the fundamental 
parameters of the accretion process such as M BH and normalized accretion rate (i.e. L/L Edd , as the 
accretion rate is not an observable), which is still not clear (see, e.g., §6 of Dietrich et al. 2002). 
There are observational reports that the ratio of the X-ray to the bolometric/optical luminosity cor- 
relates negatively with L/L Edd (Vasudevan & Fabian 2007, Kelly et al. 2008). This might account 
for the Mgll EW-£ relation as follows: as L/L Edd increases, the ratio of the X-ray to UV photons 
decreases, i.e., the heating of the Hi* region weakens and thus the Mgll EW decreases. We will 
use other emission lines such as Fe II to explore this issue in a following paper (Dong et al. 2009). 

Alternatively, we propose a physical model involving cloud properties (effectively the cover- 
ing factor). The details of the model are to be presented in the following paper (Dong et al. 2009); 
here we only present a brief description. As described in §1, Mgll line originates only from a 
thin transition layer in the partially ionized H I* region of the ionization-bounded (high-/V H ) clouds 
(see Collin-Souffrin et al. [1986] for the details). On the other hand, because low-/V H clouds are 
not massive enough to balance the radiation pressure force that is boosted due to photoelectric 
absorption as large as by about one or more orders of magnitude, they are blown out of the AGN 
BLR even at small L/L Edd . (see Marconi et al. 2008, 2009; also Fabian et al. 2006, but cf. Netzer 
2009). According to the photoionization calculation of Fabian et al. (2006, see their Fig. 1), for 
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dust-free clouds of Nu > 10 21 cm" 2 with photoionization parameter U < 1, the lower limit of the 
column density of the clouds that can survive in the BLR is approximately proportional to L/L Edd , 
as N H > 10 23 £ cm" 2 ; for dusty clouds, iV H > 5 x 10 23 ^ cm" 2 . That is, the higher the L/L Edd of an 
AGN is, the larger is the fraction of high-/V H clouds that are accumulated in the BLR. So, for Mg II 
line, it is the ionization-bounded clouds gravitationally bound in the BLR that emit it predomi- 
nantly. Thus, in the AGN ensemble, as L/L E dd increases, the number of the clouds that emit Mg II 
efficiently, in other words, the effective covering factor, decreases, and so does the Mg II EW. Thus 
the model explains the EW-£ relation naturally. The distinguishing of the two models requires 
detailed photoionization modeling, which is beyond the scope of this Letter. 
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Fig. 1 . — Plots of the equivalent width of Mg II A2800 emission line of the 2092 type-1 AGNs in 
our sample versus their Mgll FWHM, ALa(3000A), black hole mass (M B hX and Eddington ratio 
(£ = L bol /L Edd ). In panel b we also plot the apparent relation of the Mg II Baldwin effect converted 
from the best-fit logL M gii-logL 3 ooo relation with a slope of 0.91 (see Table 2), which is consistent 
with that of Dietrich et al. (2002). In panel d we also plot the best-fit linear relation in log-log 
scale between Mg II equivalent width and L/L Edd (Eqn. 1). 



-12- 




Fig. 2. — Mgll luminosity versus continuum luminosity ALa(3000A) of the type-1 AGNs in 
our sample. Objects with L/L Edd > 0.8 are denoted in blue and those with L/L Edd < 0.2 in red. 
Also plotted are the best-fit proportion relations for the full sample (green line), the subsample of 
L/L Edd > 0.8 (blue line) and the subsample of L/L Edd < 0.2 (red line), see §3.2 and Table 2. 
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Table 1 . Results of Spearman correlations and partial correlations of Mg II EW a 





FWHM 


L 






(£;FWHM) * 




' (£; m bh ) ' 


' (FWHM;£) * 


(L 30 oo;^) * 


(M m ;£) b 




0.55 


-0.18 


0.38 


-0.59 


-0.23 


-0.57 


-0.48 


0.05 


-0.02 


0.00 


Pnull C 


<le-38 


6e-16 


<le-38 


<le-38 


<le-38 


<le-38 


<le-38 


0.03 


0.3 


0.8 




0.49 


-0.22 


0.27 


-0.55 


-0.28 


-0.51 


-0.49 


0.07 


0.02 


0.04 


Pnull d 


<le-38 


8e-25 


<le-38 


<le-38 


<le-38 


<le-38 


<le-38 


0.001 


0.5 


0.07 



a For every (partial) correlation, we list the Spearman rank correlation coefficient (r s ) and the probability of the null hypothesis ( P nu n ) for the 
2092 objects in the sample. 

b (Y;Z) denotes the partial correlation between the Mgll EW and Y, controlling for Z. 

c Top panel: FWHM is measured from the model Mgll broad line; L is AL A (3000A) as measured from the power-law model of the local 
AGN continuum; The black hole masses, Mbh, are calculated using the formalism presented in McLure & Dunlop (2004); Eddington ratios 
(I = L/Lndd) are calculated assuming that the bolometric luminosity Lbol ~ 5.9AL A (3000A). 

d Bottom panel: FWHM is measured from the model H/3 broad line; L is AL A (5100A) as measured from the power-law model of the local 
AGN continuum; The black hole masses, Mbh, are calculated using the formalism presented in Collin et al. (2006, their Eqn. 7); Eddington 
ratios are calculated assuming that the bolometric luminosity L bol 9AL A (5100A). 



Table 2. Fit Parameters of log L Mg ii = HogL 3 ooo + b 



Sample 


N a 


k 


b 


Full 


2092 


0.91 ±0.01 b 


2.05 ±0.56 


Full 


2092 


1 c 


-2.03 ±0.00 


L/Lem > 0.8 


433 


0.99 ±0.02 


-1.94±0.85 


L/L Edd > 0.8 


433 


1 c 


-2.19 ±0.00 


0.2 < L/L Edd < 0.8 


1203 


0.99 ±0.01 


-1.52±0.67 


0.2 < L/L Edd < 0.8 


1203 


1 c 


-2.01 ±0.00 


L/L Edd < 0.2 


456 


1 c 


-1.83 ±0.00 


0.1<L/L Edd <0.2 


299 


0.94 ±0.04 


0.55 ±1.70 


0.1<L/L Edd <0.2 


299 


1 c 


-1.84 ±0.00 



a /V denotes the number of data points in every sample. 

b For the full sample, k is biased significantly by the se- 
lection effect, see §3.2. 

c In these fittings k is fixed to be 1 . 



